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摘要 空间 碎片 多 源 数据 融合 定 轨 软件 SPODFMD (Space debris Precise Orbit Determination Fusing Multi- 
source Data) 是 由 紫金 山 天 文 台 卫星 精密 定 轨 团 组 自主 研发 的 一 款 数值 法 编目 定 轨 软件 , 目的 是 解决 国内 空间 
: 碎片 监测 网 获取 的 多 源 数据 融合 问题 . 软件 设计 初期 , 充分 调研 了 国内 主流 空间 碎片 监测 设备 . 考虑 到 空间 碎片 
轨道 和 物理 特征 的 复杂 性 、 大 量 碎 片 轨道 更 新 对 计算 效率 的 要 求 以 及 工程 应 用 对 稳健 性 的 要 求 等 现实 因素 , 在 
定 轨 算 法 中 集成 了 具有 自主 知识 产权 的 地 球 重力 势 及 其 俩 导数 无 奇 点 快速 算法 、 无 奇 点 DTM94 (Drag Tem- 
perature Model) X ^UE] 度 模型 及 其 偏 导数 的 解析 表达 方法 、 大 偏心 率 轨道 密集 星 历 精密 快速 计算 算法 以 及 稳 
健 的 自 适 应 加 权 方 法 等 . 结合 软件 工程 设计 理论 , 目前 已 实现 了 14 种 主流 设备 获取 数据 的 任意 融合 , 并 取得 了 
型 场景 下 的 秒 级 定 轨 速度 , 适用 于 高 、 中 、 低 轨 和 各 种 偏心 率 特点 的 全 轨道 碎片 , 且 不 存在 计算 奇 点 和 极 性 奇 
点 , 经 初步 测试 证 实 了 SPODFMD 软 件 的 多 源 性 、 高 效率 、 通 用 性 和 稳健 性 效果 . 
关键 词 天 文 仪器 , 空间 碎片 , 多 源 数据 , 融合 定 轨 , 定 轨 精 度 , 计算 效率 
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1 引言 定 轨 是 组 网 观测 要 解决 的 问题 之 一 , 由 于 多 站 点 组 


eye E 网 改善 了 测量 几何 、 增 加 了 观测 频率 、 可 联合 高 精 
空间 碎片 的 轨道 信息 是 空间 态势 感知 的 重点 
全 全 度 资料 , 融合 定 轨 的 精度 优 于 单 站 定 轨 


二 关注 要 素 , 由 其 成 的 编目 数据 库 是 太空 资源 开发 编目 工作 中 的 轨道 计算 有 别 于 一 般 应 用 中 的 
三 利用、 太空 资产 碰撞 预警 、 太 空 垃圾 清理 等 应 用 的 轨道 计算 , 不 仅 要 求 轨道 计算 的 精度 足以 实现 对 空 
T 基础 , 其 重要 性 得 到 了 广泛 的 认可 . 近 20 yOR. È 间 碎 片 运动 态势 的 精确 掌控 , 还 要 求 有 较 高 的 轨道 
球 空间 碎片 专用 、 莱 用 和 可 用 监测 设备 的 数量 大 计算 效率 . 编目 工作 对 轨道 计算 效率 的 最 低 要 求 是 
WEK, 形成 以 美国 空间 监视 网 为 代表 的 多 个 监测 。 轨道 计算 的 处 理 时 间 不 能 大 于 自动 运行 设备 的 数 
网 , 设备 类 型 主要 为 无 线 电 雷达 、 光 学 望远镜 和 激 。 据 采 集 时 间 ， 以 避免 产生 数据 挤 压 的 瓶颈 效应 , 使 
光 雷 达 , 同一 目标 被 多 个 设备 探测 或 跟踪 的 概率 大 编目 过 程 难以 持续 . 以 编目 目标 2 万 、 每 天 进行 
大 增加 . 利用 多 站 点 、 多 类 型 设备 获取 的 资料 编目 。 次 轨道 更 新 计算 为 例 ， 单个 目标 的 平均 处 理 时 间 
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要 求 不 大 于 4.32 s, 对 编目 定 轨 的 计算 效率 要 求 很 高 效 性 : 软件 的 效率 在 定 轨 算法 和 编程 技巧 两 


高 . 当前 编目 工作 中 的 轨道 计算 广泛 采用 分 析 方 方面 得 到 了 优化 . 算法 方面 , 一 方面 采用 了 先 积 分 
法 , 例如 美国 采用 SGP4 (Simplified General Per- 再 插值 的 策略 替代 逐 点 积分 , 一 个 轨道 周期 积分 的 
turbation)/SDP4 (Simplified Deep space Pertur- 基点 数 随 偏心 率 增 大 由 25 弟 增 至 最 多 80 个 , 实际 构 
bation) 分 析 模 型 , 国内 编目 工作 通常 采用 拟 平均 根 造 的 基点 数 只 需要 覆盖 星 历 点 即 可 , 由 于 插值 效率 
数 法 等 . 虽然 分 析 方 法 较 好 地 解决 了 编目 工作 的 时 很 高 , 星 历 的 计算 时 间 约 等 于 插值 基点 的 积分 时 间 ， 
效 性 问题 , 但 随 着 观测 技术 的 发 展 , 它 的 模型 精度 相 比 逐 点 积分 密集 星 历 节 省 了 大 量 时 间 ; 另 一 方面 ， 
并 不 能 与 当前 的 观测 精度 (雷达 测 距 精度 优 于 20 m, 计算 一 次 积分 最 为 复杂 的 项 是 地 球 重力 及 其 偏 
光学 测 角 精度 优 于 10”) 相 匹配 , 从 而 不 能 达到 理想 数 的 计算 ， 其 耗 时 程度 随 阶 数 的 增加 呈 平 方 增长 ， 
的 定 轨 精 度 , 制约 了 相关 应 用 水 平 的 提升 , 因此 数 Mktris 系 数 变换 法 可 以 解决 这 一 问题 , 无 论 是 计算 
值 方 法 在 编目 工作 中 具有 潜在 的 应 用 价值 . 另 一 方 重力 势 还 是 重力 势 的 任意 阶 偏 导 数 , 都 可 转换 为 球 
面 , 现代 计算 机 运算 速度 的 不 断 提升 也 使 数值 方法 谐 级 数 的 求 和 , 球 谐 项 的 系数 可 事先 通过 变换 一 次 
应 用 于 编目 处 理 成 为 了 可 能 . 数值 法 应 用 于 编目 定 求 出 , 再 次 节省 了 计算 时 间 . 此 外 , 在 计算 大 气 密 度 
轨 , 应 转变 应 用 于 人 造 卫 星 定 轨 时 “一 星 一 方案 ”的 相对 于 地 固 天 坐标 的 偏 导 数 时 采用 了 解析 表达 式 ， 
> 处 理 思路 , 重点 解决 通用 性 、 时 效 性 差 的 问题 ， 具 相 比 有 限 差 分 求 导 需要 多 次 计算 大 气 密度 的 效率 
体 的 研究 方向 为 扩展 数据 类 型 、 提 升 计算 速度 、 全 要 高 . 编程 技巧 方面 , 软件 采用 了 高 效 的 排序 方法 、 
轨道 类 型 适用 以 及 海量 数据 带 来 的 计算 稳健 性 要 搜索 算法 、 滑 动 窗口 法 等 进一步 提高 计算 效率 . 
求 等 . 通用 性 : 软件 适用 于 低 、 中 、 高 轨 各 种 偏 
作者 在 以 下 方面 做 了 大 量 有 益 探索 ,包括 调 心率 的 全 轨道 空间 碎片 ， 偏 心率 通用 性 主要 通 
研 国 内 主流 空间 监测 设备 情况 (设备 类 型 、 安 装 方 过 非 均 匀 插 值 基 点 构造 方法 和 多 种 积分 器 来 实 
式 、 资 料 类 型 等 ), 调研 空间 碎片 数据 处 理 现 状 和 改 现 ， 非 均匀 插值 基点 构造 方法 是 针对 大 偏心 率 
造 要 求 (精度 和 速度 等 ), 针对 数值 法 轨道 改进 计算 轨道 利用 均匀 基点 插值 时 近地点 和 远地点 星 历 
耗 时 的 关键 环节 研发 兼顾 精度 和 效率 的 优化 算法 ， 精度 失衡 这 一 状况 提出 的 ,实际 上 适用 于 任何 偏 
针对 极 性 奇 点 和 数学 奇 点 进行 无 奇 点 化 处 理 ， 针 心率 轨道 ,可 与 任 一 单 步 法 积分 器 结合 使 用 ; 软 
对 多 源 数 据 先 验 精 度 失真 进行 自 适 应 加 权 方 法 研 件 提 供 了 4 种 积分 器 以 及 3 种 组 合 积分 器 以 适应 轨 
究 等 , 以 上 调研 结果 和 研究 结果 是 形成 SPODFMD 道 偏心 率 的 变化 ， 其 中 AC (Adams-Cowell) 积 分 
(Space debris Precise Orbit Determination Fusing 器 适用 于 偏心 率 小 于 0.1 的 轨道 , 而 RKF (Runge- 
Multi-source Data) 软 件 架构 和 定 轨 算 法 的 基础 . 本 Kutta-Fehlberg)、 RKNF (Runge-Kutta-Nystrom- 


am 


二 文 将 首先 介绍 软件 的 功能 和 特点 , 再 介绍 软件 集成 ” ”Fehlberg) 和 Everhart 适 用 于 任何 偏心 率 轨道 . 轨道 

- 的 关键 技术 , 最 后 利用 多 源 数据 进行 测试 并 根据 测 高 度 的 通用 性 主要 通过 定 轨 流 程 实现 ， 除 6 个 轨道 

试 结果 总 结 软件 的 性 能 . 根 数 外 ,高 轨 目 标 同时 改进 与 光 压 相关 的 物理 量 ， 

低 轨 目标 同时 改进 与 大 气相 关 的 物理 量 , 中 轨 目 标 

2 ”软件 的 功能 和 特点 同时 改进 与 大 气 和 光 压 相关 的 两 个 物理 量 , 通过 对 
SPODFMD 软 件 能 满足 绝 大 部 分 主流 空间 监 。 ”空间 碎片 轨道 高 度 的 区 分 实现 了 轨道 高 度 全 覆盖 . 

测 设备 获取 数据 的 融合 需求 , 并 且 具 备 多 种 特点 . 稳健 性 : 软件 的 稳健 性 主要 体现 在 无 奇 点 上 . 

Satin 软件 集成 了 两 种 地 球 重 力 势 及 其 偏 导 数 算法 ,分 

2.1 软件 功能 ( 见 附录 ) 别 为 Balmino 和 Métris 系 数 变换 法 , 它们 的 共同 特 

2.2 ”软件 特点 点 是 都 进行 了 无 奇 点 变换 ,计算 两 极 上 空空 间 点 

多 资料 : 软件 支持 14 类 设备 的 数据 输入 和 任意 ”的 引力 和 偏 导数 均 不 会 引起 定义 和 计算 的 奇 点 . 

融合 , 设备 类 型 匈 附录 . DTM94 (Drag Temperature Model) 大 气 密 度 模 型 
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也 有 具备 无 奇 点 变换 的 条 件 , 经 过 与 地 球 重力 势 相 同 
的 处 理 后 也 达到 了 消除 奇 点 的 效果 , 进一步 提升 了 
软件 的 稳健 性 . 

模块 化 : 软件 采用 模块 化 设计 ， 将 功能 独立 
有 重复 使 用 需求 的 函数 全 部 采用 模块 管理 ， 由 
此 带 来 了 程序 设计 的 灵活 性 ， 例 如 动力 学 模型 共 
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是 从 各 个 角度 同时 对 定 轨 过 程 进行 优化 , 最 终 决 定 
了 SPODFMD 软 件 的 效果 . 

3.1 大 偏心 率 轨 道 密集 星 历 精密 快速 计算 方法 
编目 工作 的 轨道 计算 普遍 涉及 密集 星 历 的 产 
生 和 计算 问题 , 星 历 点 密集 时 逐 点 积分 步 长 不 能 得 


24 个 相关 模块 ， 分 别 用 于 定 轨 / 预 报 、 低 /中 /高 
轨 、Balmino/Métris、1/2 阶 积分 器 及 以 上 4 类 配置 
选项 的 任意 组 合 , 例如 “ 定 轨 十 低 轨 十 Mktris 十 2 阶 积 
分 器 ”对 应 一 个 单独 的 模块 , 程序 的 修改 、 调 用 非 
常 便利 , 也 便于 实现 更 为 复杂 的 功能 . 

灵活 性 : 软件 的 灵活 性 是 指 多 资料 这 一 特点 在 
概念 上 的 延伸 , 软件 支持 的 14 类 设备 是 逻辑 设备 的 


到 充分 伸展 , 积分 右 函 数 计 算 次 数 大 大 增加 , 严重 
降低 数值 方法 的 计算 效率 , 避免 这 一 问题 的 有 效 技 


术 手 段 就 是 运用 内 插 方法 . 目前 产生 密集 星 历 的 多 
种 内 插 方 法 均 采 用 了 定 步 长 划分 插值 基点 的 方式 ， 
这 种 方式 适用 于 运动 特征 变化 非常 平缓 的 近 圆 轨 
道 , 对 于 编目 工作 中 经 常 涉及 的 运动 特征 变化 较为 
剧烈 的 大 偏心 率 轨 道 并 不 适用 , 易 造 成 近地点 附近 


概念 , 逻辑 设备 可 与 实际 设备 一 一 对 应 , 也 可 与 实 
际 设备 完全 不 同 , 应 根据 数据 特点 和 融合 需求 灵 沪 
设置 设备 的 类 型 . 例如 一 台 相 控 阵 雷达 ,由 于 测 距 
精度 远 高 于 测 角 精度 , 若 仅 使 用 测 距 资料 与 其 他 高 
精度 测 角 设 备 获 取 的 数据 融合 , 调用 软件 时 仅 需要 
将 相 控 阵 雷 达 的 设备 类 型 设置 为 激光 雷达 . 类 似 的 
场景 还 有 将 一 台 激 光 测 距 和 测 角 设备 作为 单 基 雷 
达 三 元 素 设备 进行 输入 等 , 采用 逻辑 设备 的 理念 使 
得 数据 融合 场景 更 加 丰富 、 人 性 化 . 
智能 化 : 智能 化 体现 在 软件 对 物理 参数 改进 方 
面 的 处 理 , 除 6 个 轨道 根 数 外 , 定 轨 时 还 尝试 同时 改 
进 与 目标 相关 的 物理 参数 , 如 果 迭 代 发 散 或 改进 值 
不 合理 则 自动 固定 其 物理 参数 重新 改进 , 整个 过 程 
无 需 人 工 干 预 . 

言 息 全 反馈 : 软件 处 理 的 所 有 环节 都 会 反馈 必 
要 的 信息 , 通过 17 类 反馈 信息 的 值 可 追溯 不 合理 的 
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输入 、 积 分 异常 、 堆 栈 异 常 等 问题 , 符合 工程 应 用 
对 软件 的 要 求 . 


3 ”关键 技术 


SPODFMD 软 件 集成 的 关键 技术 主要 有 地 球 
重力 势 及 其 偏 导数 无 奇 点 快速 计算 方法 、 无 奇 
点 DTM94 大 气 密度 模型 及 其 偏 导数 的 解析 表达 方 
法 、 大 偏心 率 轨 道 密集 星 历 精密 快速 计算 方法 以 
及 稳健 的 自 适应 加 权 方 法 等 . 这 些 算法 虽然 是 为 解 
决 某 一 问题 而 提出 , 但 在 开发 过 程 中 都 兼顾 了 定 轨 


— 


星 历 因 插 值 基点 稀疏 导致 插值 精度 低 , 而 远地点 附 
近 插值 基点 密集 而 精度 过 高 的 失衡 状况 . 
0 插值 基点 的 数量 来 改善 这 一 失衡 状况 ， 
日 相应 需 增加 数值 积分 时 间 , 这 又 易 降低 计算 效率 . 
为 兼顾 大 偏心 率 轨道 整体 的 星 历 计 算 精 度 和 计算 
效率 , 提出 一 种 非 均 匀 插 值 基 点 构造 方 法 9 引 , 具体 
为 从 已 知 星 历 的 初始 基点 开始 递 次 构造 插值 基点 ， 
待 构造 的 基点 由 上 一 个 基点 的 星 历 计算 得 出 , 基点 
的 星 历 由 数值 积分 得 出 , 基点 构造 与 数值 积分 计算 
AOS EAT, 直到 将 涵盖 所 有 插值 点 的 基点 构造 完毕 . 
非 均 匀 插 值 基点 的 构造 方法 见 (1) 式 , 基点 间隔 是 空 
间 目 标 地 心 距 的 函数 , 该 方法 可 以 与 一 般 的 单 步 积 
分 法 在 星 历 计 算 过 程 中 协调 使 用 ， 当 与 Hermite 插 
值 多 项 式 结合 使 用 时 总 效率 最 高 . 


TY 万 
An-B(u) n 


(1) 


其 中 ,和 A 了 为 第 /个 插值 基点 与 待 构造 的 第 ! + 1 个 基 
点 的 间隔 ; a 为 轨道 半 长 径 ; y 为 时 间 变 换 参 数 ; BN 
轨道 周期 变换 因子 ; Ri 为 利用 第 /个 基点 星 历 计 算 
得 出 的 地 心 距 ; 及 为 以 时 间 变 量 t 计 算 的 轨道 运动 
周期 ; 入 为 一 个 轨道 周期 内 插值 基点 的 数量 . 此 外 ， 
Y 和 为 插值 基点 构造 之 前 预先 确定 的 常量 , y 由 数 
值 实验 结果 确定 的 最 佳 普 适 数值 为 0.3, B 由 轨道 偏 
心率 e 和 Y 的 取 值 决定 , 在 基点 构造 过 程 中 也 是 一 个 


对 效率 、 精 度 、 通 用 和 稳健 的 要 求 , 它们 的 集成 就 
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3.2 “地球 重 力 势 及 其 偏 导数 无 奇 点 快速 计算 方法 

um uso 这 种 
: (1) 位 于 
两 极 上 空 Ru MU ER QURE 32 
FH XT F Hh [8] XE AVA os f e BC EY o RE HL ER BE AB SA, 
在 两 极 上 空 时 引起 分 母 为 零 . 从 物理 角度 考虑 ， 
地 球 外 部 重力 势 连续 分 布 不 存在 局 部 奇 点 ， 上述 
奇 点 非 本 性 奇 点 而 是 数学 奇 点 ， 通 过 引入 适当 的 
数学 方法 可 以 消除 . Pines 内 首先 引入 了 地 球 重 力 
势 的 一 种 无 奇 点 表达 形式 ， 将 由 勒 让 德 多 项 式 表 
达 的 地 球 重 力 势 球 谐 函 数 转换 为 由 亥 姆 霍 效 多 项 
式 表 达 ， 并 引入 地 心 距 和 方向 余 弱 取代 球 坐 标 ; 
Balmino 等 回 对 Pines 欠 的 结果 作 了 进一步 的 整理 和 
完善 , 他 们 的 基本 思想 都 是 利用 复合 求 导 法 则 将 重 
力 势 相对 于 地 固 坐 标的 一 、 二 阶 偏 导 数 转化 为 相 
对 于 地 心 距 和 方向 余弦 的 偏 导数 ,经 转化 后 的 偏 
导数 不 会 在 分 母 上 产生 纬度 余弦 的 因子 , 因此 是 
无 奇 点 的 . Hotinet 发 现任 何 球 谐 级 数 相对 于 和 矩 坐 
标的 偏 导 数 还 可 以 表示 为 一 个 球 谐 级 数 ， 于 是 将 
地 球 重 力 势 相对 于 地 固 和 珑 坐标 的 偏 导数 计算 转变 
为 一 种 系数 变换 ，Bettadpurl"1 做 了 进一步 的 研究 
和 改进 . Métris&lS mu Y Bettadpur fIPinesl 77 
法 的 结论 , 对 重力 势 展开 式 中 所 有 单个 立体 球 谐 项 
的 Cunninghaml9l 结 果 进 行 求 和 ， 并 对 和 式 中 的 立 
体 谐 函数 按 度 和 阶 进行 归并 和 整理 , 得 到 了 重力 势 
相对 于 地 固 矩 坐标 偏 导数 的 无 奇 点 紧 致 公式 , 偏 导 
o AE EE 
数 和 部 分 变量 与 原 系数 和 原 变量 之 间 存 在 变换 关 
S: e sic 
Cam > CEE”, Sim > SIO URL > LHRH, 
上 述 变量 的 具体 定义 及 变换 方法 见 下 文 (2) (3) 式 . 

SPODFMD 集 成 了 Balmino 和 Mgktris 两 种 计算 


方法 . 软件 集成 的 Balmino 方 法 对 原 公 式 中 分 开 
处 理 的 地 球 中 心 引 力 项 、 带 谐 项 和 田 谐 项 作 了 
统一 处 理 ， 表 达 简 洁 且 便于 计算 机 程式 设计 ， 其 
不 足 之 处 在 于 需要 递 推 计算 辫 姆 霍 北 多 项 式 及 其 

一 、 二 阶 偏 导 数 , 影响 计算 效率 . Mktris 方 法 实际 


为 Hotine 算 法 的 一 种 推广 , 是 迄今 为 止 能 够 实现 地 
球 重 力 势 二 阶 以 上 高 阶 偏 导数 的 唯一 方法 , 无 论 是 
计算 重力 势 还 是 重力 势 的 任意 阶 偏 导数 ， 其 计算 
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形式 是 统一 的 , 均 为 球 谐 级 数 的 求 和 ,只 是 球 谐 项 


的 系数 不 同 . 鉴于 这 些 系数 可 事先 通过 对 应 的 变 
换 求 得 , 大 大 节省 了 计算 时 间 ， 当 不 考虑 地 球 重 力 
势 谐 系数 的 变化 时 (如 忽略 潮汐 的 影响 ), 至 少 可 节 
省 30% 的 计算 时 间 [29. 

设 0U 为 重力 势 或 其 
数 , 且 有 形式 : 


-SE 


n=N1 m=0 


ef Sts E AE AVA os BA] Ek Dr fd S 


N >N L20. 


原形 式 中 , /为 地 心 引 力 常 数 ; n 和 m 分 别 为 球 
谐 级 数 展开 式 所 截取 的 阶 和 度 ， Ni 和 az 分别 为 mn 的 
下 限 和 上 限 ; 二 为 > 坐标 偏 导 的 阶 次 ,原形 式 中 其 值 
为 零 ; Cam 和 Sum 为 归 一 化 的 球 谐 系 数 ， 刀 sm 为 归 
一 化 的 赫 姆 霍 效 多 项 式 ， HA ym = cos" op cos mA, 
om = cos" yp sin mA, 入 和 2 分 别 为 目标 点 的 经 度 和 
纬度 . MIU TEDSE-T XB A B (a, y, 2) RAE (i, j, K) TI 
偏 导 数 公 式 为 : 


(2) 


= 


Oititk ul m 
oot) 
Ox'0ylOz* r > r A 

?一 Ni m=0 
(ee + SG Ji Pa.) Hym, (3) 
中 ,py = uu. L = Ltk, Ni = Niti+j+k 


Ng — Nz i jk, Con” 和 Ss 的 变换 公式 及 


符号 的 含义 参见 文献 [9]. 


3.3 ” 自 适 应 加 权 方 法 

\ 有 不 同 测量 精度 的 多 类 型 资料 融合 时 需 重 
点 考虑 权 的 分 配 问 题 , 常见 的 做 法 是 根据 资料 的 先 
验 精度 加 权 , 先 验 精度 失真 易 造成 轨道 改进 质量 的 
PRE. 考虑 到 精度 失真 是 一 种 常见 现象 , 提出 了 一 
种 稳健 的 自 适应 加 权 方 法 , 自 适 应 权 是 在 固定 权 轨 
道 改 进 收 敛 的 基础 上 利用 资料 残 差 重 新 分 配 权重 
直至 收敛 的 过 程 , 其 特点 在 于 : (1) 自 适应 加 权 定 轨 
利用 了 固定 权 轨 道 改进 收敛 的 结果 作为 初 根 数 , 因 
此 启动 前 会 默认 先 执行 固定 权 定 轨 , 这 样 处 理 可 避 
免 直接 使 用 自 适 应 权时 ( 自 适 应 权 由 初 轨 预报 值 与 
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资料 值 的 残 差 确定 ) 过 度 依赖 初 根 数 的 精度 , 这 种 ”4 ”测试 与 分 析 
依赖 性 使 得 轨道 改进 偏离 真 值 的 概率 大 大 增加 , NEA B 
多 源 观测 数据 有 两 层 含义 , 一 是 指 来 源 于 多 个 
noo E RU, 二 是 指 多 种 数 据 类 型 ,两 者 满足 其 一 即 可 称 
观测 资料 , 利用 固定 权 轴 道 改进 最 后 一 次 迁 代 刚 除 ”为 融合 定 轨 ， 以 下 将 通过 算 例 来 测试 软件 的 性 能 
野 值 后 的 全 部 资料 参与 自 适 应 定 轨 , 这 样 处 理 既 对 ae os ec a MERE UI S 
精度 失真 资料 的 权 进 行 了 一 定 程度 的 纠偏 , 又 保证 
了 自 适应 迭代 的 稳健 性 , 通常 自 适应 权 定 轨 仅 需要 41 ” 单 站 资料 定 轨 
少量 的 迭代 即 可 收敛 , 不 会 显著 增加 计算 时 间 . 单 站 资料 指 一 台 “ 逻 辑 设 备 2> 同 时 获取 的 测 轨 
无 亲 点 DiTML94 大 气 密度 病 型 及 其 仿 导 数 的 “资料 , 例如 雷达 可 同时 获取 测 角 和 测 距 三 元 素 资料 , 


解析 表达 方法 

在 低 轨 目标 编目 定 轨 中 需要 考虑 大 气 阻 力 摄 
动 ， 对 于 轨道 特别 低 或 者 面 质 比特 别 大 的 低 轨 目 
bs, 大 气 阻力 摄 动 已 与 地 球 主 要 带 谐 项 .项 摄 动 
相当 , 因此 为 了 保证 定 轨 迭 代 过 程 的 收敛 , 还 必须 


有 的 雷达 还 能 获取 包括 测速 在 内 的 四 元 素 资料 . 同 
时 获取 的 三 元 素 或 者 四 元 素 都 是 单 站 资料 , 大 部 分 
单 站 资料 都 是 由 单 台 实际 设备 获取 的 , 也 有 部 分 单 
站 资料 涉及 到 收发 分 置 的 两 台 实 际 设备 , 例如 双 基 
雷达 的 距离 和 资料 以 及 无 线 电 时 差 和 频 差 资料 , 此 
时 一 台 “ 逻 辑 设 备 " 由 两 台 实 际 设备 构成 . 包含 同类 


计算 大 气 阻 力 相 对 于 目标 位 置 和 速度 矢量 的 
数 . SPODFMD 采 用 了 DTM94 热 大 气 模 型 41， 其 
为 三 维 模型 具有 精度 高 的 优点 , 但 由 于 也 采用 球 
坐标 作为 基本 变量 , 因此 模型 本 身 及 其 相对 于 地 固 
和 矩 坐标 的 偏 导 数 也 存在 极 性 奇 点 ; 男 一 方面 , 目前 
除 Harris-Priester 大 气 模 型 已 有 成 熟 的 偏 导 数 解析 
表达 式 外 ， 其 他 大 气 模型 均 因 过 于 复杂 难以 求 导 
利用 差分 计算 偏 导数 是 常用 的 蔡 代 方 法 , 然而 当 模 
型 不 连续 或 不 可 微 时 , 差分 计算 结果 失真 . 

鉴于 DTM94 原 模型 中 存在 勒 让 德 多 项 式 , 同 
样 借鉴 Pines 引 消除 地 球 重力 势 球 谐 表 达 式 奇 点 所 
的 方法 , 用 目标 地 心 距 以 及 方向 余弦 取代 球 坐 
标 作为 DTM94 大 气 模型 的 基本 变量 , 变换 后 的 大 
气 密度 模型 是 无 奇 点 的 ， 且 使 得 寻求 大 气 密度 相 
对 于 空间 目标 地 固 矩 坐标 的 解析 表达 成 为 可 能 
借鉴 Pines 欠 提出 的 复合 求 导 法 则 ， 提 出 了 大 气 密 


局 导 
其 


P 
— 


RE PAART FHA SURE AE RI RA RI: 
ap 109 | (8p 1 记 ap \ 
OX; rð E 7 x^) 4» wi 


其 中 ,了 = 1,2,3, p = p(z,y,z)?3yx. y. 2z 的 连续 
可 微 函数 ， HAX: = T, X2 = 4; X3 = 2, Ts 01. 
05. 032) WA H ES ESL HR CHEER S] A FX Y, 
2Z 轴 上 的 方向 余弦 . 相 比 差分 求 取 偏 导数 的 方法 , 解 
析 表 达 式 具有 准确 、 稳 健 和 高 效 的 优势 
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型 数据 的 单 站 资料 , 例如 激光 测 距 、 光 电 测 角 这 类 
资料 由 于 仅 包含 一 种 元 素 或 者 精度 相同 的 两 种 元 
A, 定 轨 时 等 权 分 配 权 重 即 可 , 而 三 元 素 和 四 元 素 
资料 包含 了 不 同 精度 的 测 角 、 测 距 和 测速 3 类 数据 ， 
即便 来 源 于 同一 设备 , 也 涉及 不 同 数 据 类 型 的 融合 
问题 , 空间 碎片 多 源 观测 数据 融合 定 轨 软 件 采用 马 
尔 可 夫 估 计 理 论 , 利用 数据 的 先 验 精度 确定 各 元 素 
的 权重 , 解 算 结果 稳定 准确 . 以 下 通过 例子 来 说 明 
单 站 资料 的 定 轨 能 力 , 各 算 例 的 数据 类 型 、 时 间 跨 
度 和 数据 量 等 关键 信息 列 于 表 1, 利用 这 些 资料 的 
定 轨 性 能 主要 指标 列 于 表 2,， 采 用 的 主要 动力 学 模 
型 和 公共 配置 参数 见 表 3. 
算 例 1-5 是 较为 典型 的 雷达 获取 空间 碎片 的 数 
据 条 件 ， 数 据 的 时 间 跨 度 约 2-5 d， 每 天 能 探测 1- 
2 次 , 整体 来 说 符合 空间 碎片 稀疏 性 的 特点 . PAC 
据 条 件 下 轨道 改进 收敛 通常 需要 较 多 的 迭代 次 数 ， 
由 表 2 的 统计 结果 可 知 以 上 算 例 在 SPODFMD 软 件 
中 的 迭代 次 数 为 3-8 次 , AC、Everhart 和 RKF 积 分 
器 的 平均 处 理 时间 分 别 为 0.91 s、1.62 s 和 1.78 s, 积 
分 器 的 计算 效率 以 AC 10 最 佳 , Everhart 11 和 RKF 
8 (7) 弟 次 降低 . AC 积 分 器 的 高 效 众 所 周知 ， 然 而 
可 变 步 长 的 单 步 法 积分 器 Everhart 和 RKE 的 计算 
效率 接近 AC 积 分 器 也 是 因为 采用 了 先 积 分 再 插值 
的 策略 , 而 星 历 精度 的 保障 则 主要 归功 于 关键 技术 
一 一 非 均 匀 插 值 基 点 构造 方法 的 应 用 . 该 方法 可 与 


= 
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任何 一 种 单 步 法 积分 器 结合 使 


D 


心率 轨 


道 星 历 精度 和 计 和 


1 


1, 达到 兼顾 各 
效率 的 目的 , 这 一 效果 不 


fimi 


H 


论 是 定 步 长 积分 器 AC 还 是 变 步 长 积分 器 逐 点 
均 不 能 实现 . 


表 1 单 设备 (SE) 获 取 的 数据 资料 
Table 1 Data from Single Equipment (SE) 


Case No. Object Orbit Orbit Equipment Data Data time Data 
ID type eccentricity type type span/h volume 
SE-1 1646 LEO 0.02057 GBR A,E,d 31.772 481 
SE-2 16908 LEO 0.00114 GBR A,E,d 50.316 534 
SE-3 16908 LEO 0.00114 GBR A,E,d 50.357 803 
SE-4 37210 GEO 0.00050 GBR A,E,d,d 92.842 632 
SE-5 38953 GEO 0.00058 GBR A, E, d, d 92.836 632 


Notes: A, E, d, and d represent azimuth, elevation, range, and radial velocity respectively. 
GBR is the abbreviation of Ground-Based Radar. 


表 2 SE 定 轨 算 例 性 能 


Table 2 Performance of SE orbit determination cases 


Case Computation time/s Accuracy Iteration Precise 

Data element . . 

No. RKF 8 (7) Everhart 11 AC 10 /m times orbit 
A, E,d 1.45 1.23 0.53 79 6 

SE-1 A, E 2.07 1.86 1.08 2059 7 CPF 
d 2.06 1.84 1.03 79 7 
A,E,d 1.95 1.67 0.61 144 6 

SE-2 A,E 2.23 2.03 1.03 6272 5 CPF 
d 2.45 2.20 1.17 186 6 
A,E,d 1.76 1.53 0.56 22 5 

SE-3 A,E 1.79 1.64 0.97 3610 3 CPF 
d 2.06 1.87 1.03 23 4 
A,E,d,d 1.31 1.29 0.97 7265 6 

SE-4 SP3 
A,E,d 1.19 1.08 0.89 18531 4 
A,E,d,d 1.22 1.23 1.08 5234 8 

SE-5 SP3 
A,E,d 1.72 1.70 1.01 25970 5 


Notes: CPF (Consolidated Prediction Format) and SP3 (Standard Product #3) are two kinds of 
precise orbits released by ILRS (International Laser Ranging Service) and IGS (International GNSS 
Service) respectively; The accuracy of each orbit determination case is assessed by the RMS (Root 


Mean Square) of the position differences between computation orbit and precise orbit within one 


period centered by the orbit determination epoch. 
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数据 条 件 是 影响 定 轨 精度 的 主要 因素 ,利用 观测 资料 在 单 站 轨道 改进 中 发 挥 了 决定 性 作用 , 在 
同一 时 刻 获 取 多 维 资料 中 的 不 同 元 素 及 其 组 合 各 元 素 先 验 精度 较为 准确 的 ' 
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RI SPB 定 轨 的 力学 模型 及 公共 参数 


Table 3 Force models and public parameters of SE orbit determination 


Item Model and recommended value 
GRIM4.S4; 20 x 20 for LEO objects 
and 8 x 8 for GEO objects 


Degree and order for gravity vector 


Degree and order for gravity gradient tensor 4x4 
Geopotential derivatives algorithm Métris 
Nutation IAU 1980, 4 terms 
Atmospheric density and derivatives D'TM94 and analytic formula 
Lunisolar attraction Jean Meeus analytic formula 
Weighting method Fixed weight 


一 


青 况 下 , 融合 多 维 资料 


定 轨 ， 单 站 算 例 1-3 仅 利用 了 雷达 测 距 元 素 的 轨 的 轨道 改进 效果 一 般 不 进 于 单 维 资 料 - 
道 改进 收敛 且 定 轨 精 度 达 到 了 百 米 甚 至 以 下 量 42 多 源 数据 定 轨 


2k, MAI F 


了 雷达 测 角 元 素 的 轨道 改进 精度 在 公 


多 源 数据 融合 要 求 同 一 时 间 段 内 能 获取 到 来 


里 级 别 ， 同 


NRI RUMER-UURISSGEUUR — 自 于 不 同 设备 的 观测 数据 , 当 组 网 观测 或 者 协同 观 


精度 不 低 于 仅 利 用 高 精度 测 距 元 素 的 定 轨 精度 ; 测 时 会 产生 多 源 数 据 融 合 定 轨 的 需求 . 以 下 用 典型 


单 站 算 例 4-5 针 对 高 轨 目 标的 雷达 观测 增加 了 测 。 多 站 场景 来 测试 软件 的 性 能 和 融合 效果 , 表 4 列 出 


E 


Bi 


速 元 素 ， 仅 利用 测 角 或 测 距 元 素 的 定 轨 误 差分 别 了 7 个 融合 场景 的 数据 特点 , 表 5 为 对 应 场景 的 定 轨 
为 25 km 和 18 km, 增加 高 精度 测速 元 素 的 融合 定 轨 性 能 关键 指标 ， 表 6 为 主要 的 动力 学 模型 及 公共 配 
误差 分 别 降 至 7 km 和 5 km. 以 上 算 例 体现 了 高 精度 置 参 数 . 


表 4 多 设备 (ME) 获 取 的 数据 资料 
Table4 Data from Multi Equipment (ME) 


Case No. Object Orbit Orbit Equipment Data Data time Data 
ID type eccentricity type type span/h volume 

ME-1 37210 GEO 0.00050 1SBO & 3GBO a, 0 16.230 1258 
ME-2 38001 GEO 0.00094 1SBO & 1GBO oo 0 48.716 88 
ME-3 16908 LEO 0.00114 2GBR A,E,d 69.392 1337 
ME-4 7646 LEO 0.02057 2GBR A, E,d 52.505 1498 
ME-5 46469 LEO 0.00061 3GBL & IGBO A, E,d 111.013 2743 
ME-6 1328 LEO 0.02580 2GBR A, E,d 52.218 1992 
ME-7 7647 LEO 0.02182 3GBR A, E,d 69.379 1014 


Notes: SBO, GBL, and GBO are abbreviations of space-based optical telescope, ground- 
based laser, and ground-based optical telescope respectively. Therefore, 1GBO represents one 


ground-based optical telescope. 
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#5 ME 定 轨 算 例 性 能 


Table 5 Performance of ME orbit determination cases 


Case No. Computation time/s Accuracy/m Iteration times Precise orbit 


ME-1 0.80 
ME-2 0.73 
ME-3 1.93 
ME-4 1.54 
ME-5 1.30 
ME-6 1.70 
ME-7 1.39 


433 2 EPH 
4748 7 SP3 
36 6 CPF 
40 8 CPF 
49 4 CPF 
26 3 CPF 
23 2 CPF 


Notes: EPH is a type of GNSS data that stores precise ephemerides of a satellite. 


#6 ME 定 轨 的 力学 模型 及 公共 参数 


Table 6 Force models and public parameters of ME orbit determination 


Item 


Model and recommended value 


Degree and order for gravity vector 


GRIM4_S4, 20 x 20 for LEO 
objects and 8 x 8 for GEO objects 


Degree and order for gravity gradient tensor 4x4 


Geopotential derivatives algorithm 


Nutation 


Atmospheric density and derivatives 


Lunisolar attraction 
Integrator 


Weighting method 


Métris 
IAU 1980, 4 terms 
DTM94 
Jean Meeus analytic formula 
Adams-Cowell (10) 


Fixed weight 


融合 算 例 1-7 展 示 了 多 源 数 据 融合 定 轨 的 效 


GE E 


基 


Er 
Lr 
RÀ 


=% 


, 表 5 显 示 高 轨 目 标 37210 利 用 1 台 天 基 光 学 和 3 人 台 


光学 融合 的 精度 达到 了 433 m， 相 比 利 用 单 


雷达 四 元 素 资 料 7 km 的 定 轨 精 度 提高 了 一 个 


以 上 ; 高 轨 目 标 38091 利 用 1 台 天 基 和 1 台地 


学 的 融合 精度 优 于 5 km; 低 轨 目标 16908、 


7646 和 46469、1328 和 7647 的 融合 精度 均 优 于 50 m. 


ap" 
虽然 


HB 


依据 以 上 少量 算 例 不 能 直接 得 出 多 站 融合 精 


及 局 
普遍 
引起 
7 利 月 


于 单 站 的 结论 , 单 从 数值 上 分 析 , 融合 的 精度 
高 于 单 站 , 这 个 提升 主要 由 数据 几何 条 件 改善 
, 数据 质量 和 数量 是 弱 相 关 因 素 . 融合 算 例 1- 


站 融 


合 定 轨 处 理 的 数据 量 相 比 单 站 有 所 增加 , 但 更 


AC 积分 器 的 平均 计算 时 间 为 1.34 s， 虽 然 多 


加 丰富 的 数据 使 得 迭代 次 数 减 少 , 总 体 而 言 计算 时 
间 仅 有 少许 增加 , 仍 处 于 秒 级 水 平 . 


4.3 ”动力 学 模型 测试 

SPODFMD 软 件 支 持 70 阶 以 内 地 球 重 力 势 及 
其 偏 导 数 的 计算 , 并 可 扩充 附加 潮汐 、 地 球 反 照 辐 
射 压 、 行 星 引 力 等 摄 动力 . 针对 低 轨 目标 7647 和 高 
轨 目 标 37210， 图 1 和 图 2 分 别 展 示 它 们 利用 简约 动 
力学 模型 和 精密 动力 学 模型 进行 轨道 确定 并 预报 
轨道 的 位 置 误差 (为 表述 方便 ， 由 简约 动力 学 模型 
定 轨 及 预报 的 轨道 称 为 F 轨 道 ,由 精密 动力 学 模型 
定 轨 及 预报 的 轨道 称 为 P 轨 道 ), 在 定 轨 期 内 7647 目 
标的 F 轨 道 与 P 轨 道 基 本 重合 , 预报 至 2 d 两 者 的 偏 
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差 约 50 m, 预报 至 6 d 两 者 的 偏差 约 230 


标 37210 整 个 定 轨 及 预报 期 内 F 轨 道 与 P 轨 道 几 乎 


ChinaXiv 合 作 期 刊 


间 碎 片 多 源 数 据 融合 定 轨 软 件 SPODEFMD 简 介 5 期 


m; 高 轨 目 标 取 8 x 8 阶 ,地球 重力 势 偏 导数 计算 取 4 阶 ， 章 动 


模型 取 4 项 ， 考 虑 日 月 引力 , 低 轨 目标 考虑 大 气 阻 


重合 . EARN 


了 精密 动 


片 编目 定 轨 . 基于 这 


， 两 个 例子 均 显示 利用 
力学 模型 的 轨道 改进 时 间 相 比 利 用 简约 动力 学 模 


型 增加 了 约 30% 的 计算 时 间 . 己 知 编目 目标 的 轨道 
更 新 频率 约 1-2 d 1 次 , 短期 内 简约 动力 学 模型 的 
定 轨 和 预报 精度 与 精密 模型 的 差距 不 显著 , 因此 不 
推荐 牺牲 计算 效率 的 精密 动力 学 模型 用 于 空间 碎 
WU, SPODFMD 软 件 在 测 


JJ, 高 轨 目 标 考虑 太阳 光 压 , 采用 以 上 简约 模型 即 
可 取得 表 2 和 表 5 中 的 定 轨 效 果 . 图 1-2 用 于 测试 的 
精密 动力 学 模型 的 配置 如 下 : 低 轨 目标 取 地 球 引 力 
场 20x20 阶 , 中 高 轨 目 标 取 8x8 阶 , 地 球 重力 势 偏 导 
数 计 算 取 4 阶 , 章 动 模型 取 108 项 , 考虑 日 月 、 行 星 
引力 和 地 球 反 照 辐射 压 , 低 轨 目标 考虑 大 气 阻 力 和 
WY GRRE. TELA. Aa), 高 轨 目 标 考虑 太阳 


试 算 例 时 , 均 采 用 了 简约 的 动力 学 模型 
如 下 : 低 轨 目标 取 地 球 引 力 场 20 x 20 阶 


800 
| ID: 7647 
Altitude: 959 km 


JEJE. 


， 配置 方法 
, 中 高 轨 目 


600 | Data: 3GBR / 69.379h / 1014 volume /9 pass 


Position error/m 
CD A ol 
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00 L Orbit determination time: 1.39s (F) vs 1.76s (P) 


Orbit prediction/d 


Eli 低 轨 


标 分 别 利 


简约 动力 学 模型 和 精密 动力 学 模型 定 轨 并 预报 轨道 的 位 置 误差 


Fig.1 Orbital position errors of a LEO object determined and predicted by using a simplified dynamic model and a precise 


dynamic model 
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Orbit determination time: 0.80s (F) vs 1.08s(P) 
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图 2 高 轨 


标 分 别 利 


简约 动力 学 模型 和 精密 动力 学 模型 定 轨 并 预报 轨道 的 位 置 误差 


Fig.2 Orbital position errors of a GEO object determined and predicted by using a simplified dynamic model and a precise 


dynamic model 
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Introduction to Orbit Determination Software SPODFMD for 
Multi-source Data Fusing of Space Debris 


DU Jian-i^? | XU Jin? | YANG Dong? CAO Zhi-bin22 MA Jian-bo!? 


(1 Purple Mountain Observatory, Chinese Academy of Sciences, Nanjing 210023) 
(2 Key Laboratory of Space Object and Debris Observation, Chinese Academy of Sciences, Nanjing 210023) 


Asstract The data-fusion orbit determination software SPODFMD (Space debris Precise Orbit Deter- 
mination fusing Multi-source Data) is developed with the purpose of cataloging space debris, characterized 
by integrating pure numerical integrators and given birth in the lab of SPOD (Satellite Precise Orbit De- 
termination) at PMO (Purple Mountain Observatory). At the initial stage of the birth, developers fully 
investigated national space debris monitoring equipment, fully understood the orbit and physical com- 
plexity of numerous space debris, fully recognized the requirement for high efficiency when renewing orbits 
for lots of debris, and fully considered the robustness requirement in engineering. In addition, SPODFMD 
integrates several algorithms with indigenous intellectual property rights, including a rapid and non- 
singular algorithm for computing the Earth's gravitational potential and its derivatives, analytical and 
non-singular expressions for the DTM94 thermosphere model and its derivatives, an accurate and effi- 
cient method for calculating dense ephemeris of high-eccentricity orbit, and a robust adaptive weighting 
method. Integrating these advanced algorithms by applying software engineering theory, SPODFMD can 
freely fuse observation data from 14 kinds of equipment, achieving a second-level computation efficiency 
in a majority of typical orbit determination cases. Moreover, the software shows similar performance for 
GEO (Geosynchronous Equatorial Orbit), MEO (Medium Earth Orbit), LEO (Low Earth Orbit), and 
HEO (Highly Elliptical Orbit) objects without computation and pole singularity. It is tested and proved 
to be multi-source, efficient, general, and robust. 


Key words astronomical instrumentation, space debris, multi-source data, data-fusion orbit determina- 
tion, orbit determination accuracy, computation efficiency 
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附录 


表 7 SPODFMD 软 件 功 能 简 述 
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Table7 Brief introduction of the software SPODFMD's function 


Function 


Description 


Notes 


Initial orbital 


elements 


Equipment type 


Target type 


Earth gravity 


Appliable to 3 types of initial orbital elements: TLE, f, F, 


and quasi-mean orbital elements 


Appliable to 14 types of space-surveillance equipment: 
1. Radar yielding 3-dimension data [A, E, d]; 

2. Radar yielding 4-dimension data [A, E, d, d]; 

3. Bistatic radar yielding 3-dimension data [A, E, 2d]; 

4. Bistatic radar yielding 4-dimension data [A, E, 2d, d]; 
5. Alt-az photoelectric telescope yielding azimuth and 
altitude angular data [A, E]; 

6. Alt-alt photoelectric telescope yielding longitude and 
latitude angular data [L, B]; 

7. Equatorial photoelectric telescope yielding hour and 
declination angular data [t, o]; 

8. Photoelectric telescope yielding right ascension and 
declination angular data [«, 6]; 

9. laser ranging data [d]; 

10. bistatic laser ranging data [2d]; 

11. space-based telescope yielding right ascension and 
declination angular data [&, 6]; 

12. ship-based radar yielding 3-dimension data [A, E, d]; 
13. ship-based radar yielding 4-dimension data [A, E, d, d]; 
14. ship-based telescope yielding right ascension and 


declination angular data [a, 5]; 


Applicable to space objects with various eccentricities in all 


orbital regimes around the Earth 


Supporting the user respectively in setting the degree and 
order of the truncated spherical harmonic expansion of 
geopotential within a range of 4 to 70 when computing 


gravity and its partial derivatives 


Other types can be converted 
to one of the supported types 
ahead 


Equipment is a kind of logical 
concept, indicating one or more 
pieces of actual equipment that 
can obtain one of 14 types of 
data. 

A piece of conventional equip- 
ment involves one station, and 
a bistatic involves two stations. 
Equipment without marking 
out the carrier’s type indicates 


a ground-based one. 


Apart from six orbital ele- 
ments, the software also cor- 
rects the physical parameters 
associated with atmospheric 
solar 


drag or/and pressure 


when processing rich data. 


The recommended degrees and 
orders for gravity computation 
use are 20 x 20 and 8x8 for 
LEO and GEO objects re- 
spectively, whereas the rec- 
ommended degree and order 
for partial-derivatives compu- 


tation use is 4x4 for all objects. 
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表 7 续 
Table7 Continued 
Function Description Notes 
Geopotential Providing two algorithms of Balmino and The more efficient algorithm Métris applies to 
algorithm Métris for the computation of gravity and the geopotential with fixed spherical harmonic 


Nutation term 


Integrator 


Weighting 
method 


Orbit 
determination 


epoch 


Residual file 


Feedback 


its partial derivatives 


Supporting the user in setting the number 
of terms of the IAU 1980 nutation model 


Providing 7 integrating ways as well as 
optional orders of integrators: 

1. AC: 8~12 orders and PE or PCE 
method; 

. RKF: 6(5), 7(6), 8(7), and 9(8) orders 
. RKNF: 5(6), 6(7), and 7(8) orders 

. Everhart: 11 and 15 orders 
AC+RKF: set orders respectively 
AC+RKNF: set orders respectively 


Noa o&Roo nw 


. AC+Everhart: set orders respectively 


Providing two weighting methods: fixed and 


adaptive 


Supporting the user in setting the epoch of 
orbit determination within the data range 


by utilizing an adjusting factor 


Providing the switch to turn on or turn off 


printing the residual file 


Providing the function of outputting the 


feedback of orbit determination 


coefficients. 


The recommended value is 4. 


The AC integrator is suitable for objects of 
quasi-circular orbits (eccentricities below 0.1), 
being superior to other integrators in efficiency. 
The one-step integrators, such as RKF, RKNF, 
and Everhart, are applicable to all eccentricities. 
The AC integrator can be utilized along with 
a one-step integrator if the eccentricity floats 
around 0.1 in the iteration procedure. The rec- 
ommended orders are 10, 8(7), 7(8), and 11/15 
for integrators of AC, RKF, RKNF, and Ever- 


hart respectively. 


The fixed weight is determined according to each 
piece of the data’s prior precise and keeps un- 
changed, while the adaptive weight adjusts ac- 
cording to the posterior precise in each iteration 


procedure till converging. 


The adjusting factor is a value between 0 and 1, 
and the epoch of orbit determination equals to 
the factor multiplied by the data span plus the 


data’s start time. 


It is a suggestion that turning off the printing 
function to reduce time and memory consump- 
tion for daily operation, and turning it on when 
manual interpreting is required for few anoma- 


lous cases. 
The feedback includes 17 types of informa- 
tion, involving astringency tips, input parame- 


ters check, stack exception tips, etc. 
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